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We explain the mechanism leading to directed chaotic transport in Hamiltonian systems with 
spatial and temporal periodicity. We show that a mixed phase space comprising both regular and 
chaotic motion is required and derive a classical sum rule which allows to predict the chaotic trans- 
port velocity from properties of regular phase-space components. Transport in quantum Hamiltonian 
ratchets arises by the same mechanism as long as uncertainty allows to resolve the classical phase- 
space structure. We derive a quantum sum rule analogous to the classical one, based on the relation 
between quantum transport and band structure. 



Stimulated by the biological task of explaining the 
functioning of molecular motors, the study of ratchets |l]] 
has widened to a general exploration of "self-organized" 
transport, i.e., transport without external bias, in nonlin- 
ear systems Q. Along with this process, there has been a 
tendency to reduce the models under investigation from 
realistic biophysical machinery to the minimalist systems 
customary in nonlinear dynamics. External noise, for ex- 
ample, which originally served to account for the fluctuat- 
ing environment of molecular motors, has been replaced 
by deterministic chaos. This required to include inertia 
terms in the equations of motion, thus leaving the regime 
of overdamped dynamics and leading to deterministic in- 
ertia ratchets with dissipation It is then a con- 
sequent but radical step to abandon friction altogether. 
Indeed, transport in Hamiltonian ratchets was observed 
numerically if all symmetries were broken that generate 
to each trajectory a countermoving partner [@J1. 

As a parallel development, the desire to realize ratchets 
in artificial, nanostructured electronic systems, required 
to consider quantum effects [j^,[3| . Quantum Hamiltonian 
ratchets, however, have been studied only in the frame- 
work of one-band systems where no transport occurs Q . 

In this paper we explain how a Hamiltonian ratchet 
works. We rely on methods which — although well estab- 
lished in studies of deterministic dynamics — have never 
before been applied to ratchets. We derive a classical 
and an analogous quantum sum rule for transport allow- 
ing the following conclusions: (i) Directed transport is 
a property associated with individual invariant sets of 
the dynamics. A necessary condition for non-zero trans- 
port is a mixed phase space with coexisting regular and 
chaotic regions, (ii) Transport in chaotic regions can be 
described quantitatively by using topological and further 
properties of adjacent regular regions only, (iii) Quantum 
transport persists for all times and approaches the clas- 
sical transport when K is small compared to the major 
invariant sets of the classical phase space. 
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FIG. 1. (a) Spatial distribution P(x) of a continuously 
driven system || after 20, 000 time periods showing the di- 
rected transport in a Hamiltonian ratchet. Initially, 10 4 tra- 
jectories were started at random on the line p — 0, x 6 [0, 1) 
in the chaotic sea. (b) Poincare section p vs x of a unit cell 
at integer times showing the main chaotic sea, the upper and 
lower limiting KAM-tori, and the major embedded regular 
islands. (c,d) As (a,b), but for the kicked Hamiltonian (^) 
showing a much more pronounced directed transport. 

We consider a Hamiltonian of the form H(x,p,t) = 
T(p)+V(x, t), where T(p) is the kinetic energy. The force 
—V is periodic in space and time, V'(x+l,t) — V'(x, t + 
1) = V'(x, i), and has zero mean L dt f Q dx V'(x, t) = 0. 
Usually directed transport is demonstrated by following 
selected trajectories over very long times |^||| or an en- 
semble of trajectories which generates spatial distribu- 
tions as shown in Fig. [l]a,c. While this is easily imple- 
mented numerically, it gives no clue about the origin of 
the transport (but see Ref. M). Instead, we shall exploit 
the periodicity of the dynamics with respect to space and 
time and analyze transport in terms of the invariant sets 
of phase space, reduced to the spatio-temporal unit cell 
x, t £ [0,1). For any finite invariant set M we define 
ballistic transport as phase-space volume times average 
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velocity expressed as 

dt dx I dpxM{x,p,t) 
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where \m is the characteristic function of M. Transport 
is additive for the union of two or more disjunct invariant 
sets, i.e., for M = 1J- Mi, with M, n M = for all i^j, 
we have 
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This sum rule for Hamiltonian transport has far-reaching 
consequences to be discussed in the following. 

For a generic Hamiltonian system, phase space is mixed 
and comprises an infinite number of minimal invariant 
sets of different types. For the sake of definitcness we 
will restrict the following discussion to the most interest- 
ing case of a chaotic region containing embedded regu- 
lar islands (Fig. |l|b). In any of these invariant sets the 
time-averaged velocity V{ is the same for almost all initial 
conditions (assuming ergodicity for chaotic components) . 
Hence, for a chaotic region, r c h = A^ v c h with A c h de- 
noting its area in a stroboscopic Poincare section. For an 
embedded island we have Tj = Ai Vi where Ai includes 
the areas of the narrow chaotic layers inside the island 
and of the infinite hierarchy of island chains surrounding 
it because all these invariant sets share the same mean ve- 
locity. This velocity Vi is identical to the rational winding 
number w; = Xi/U of the stable fixed point at the center 
of the island. In extended phase space, this corresponds 
to a shift of the island by Xi spatial after ti time periods. 
Typically, the chaotic set is bounded from above and be- 
low by two non-contractible KAM-tori p a ,b(x, t) enclosing 
the spatial unit cell. Treating the phase-space region in 
between as the global invariant set M appearing on the 
l.h.s. of Eq. (||), its transport tm is obtained from Eq. (|l|) 
as (T% - (T) b with {T) a . b = J Q dt J Q dxT(p a . b (x,t)) de- 
noting averages of the kinetic energy over the tori. Using 
the sum rule @ we can now express transport of the 
chaotic region in terms of its adjoining regular compo- 
nents (KAM-tori and islands) as 



A^ch = (T) a - (T) b 



Ai Wi . 
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This is our main result on classical transport in Hamilto- 
nian ratchets. Not only does Eq. (Q) provide an efficient 
method to determine the chaotic drift velocity, it also ex- 
presses the simple principle generating directed ballistic 
motion: Decomposing phase space into different invariant 
sets, these will in general have average velocities different 
from each other and also different from zero but related 
by the sum rule (^) . Therefore a necessary condition for 
directed chaotic transport in Hamiltonian ratchets is a 
mixed phase space. 



Levy flights |T(| are a characteristic feature of chaotic 
motion in a generic mixed phase space and indeed they 
were observed in Hamiltonian ratchets |^||. They re- 
flect the slow exchange between subsets of a chaotic re- 
gion, separated by leaky barriers |l(|. As these subsets 
are not invariant, their contributions to Eq. (^|) are con- 
tained in the contribution of the chaotic invariant set. 
Levy flights lead to power-law tails in spatial distribu- 
tions. For example, the asymmetric shapes of the peaks 
visible in Figs. [l]a,c can be attributed to these tails. No- 
tably, in Fig. |l|c one clearly sees a mean transport to the 
right although the same data shows no indication of a 
power-law tail in this direction. We stress that the sum 
rule allows to predict the mean velocity of chaotic trajec- 
tories without any reference to such details of the chaotic 
dynamics. This suggests that Levy flights are not a nec- 
essary element of the mechanism of chaotic transport in 
Hamiltonian ratchets. 

In Ref. § it was shown that a necessary condition 
for directed transport is the breaking of all symmetries 
which to each trajectory generate a countermoving part- 
ner. For a chaotic set invariant under such a symmetry, 
this is in agreement with Eq. (||) because then the r.h.s. 
vanishes identically. However, chaotic sets can also oc- 
cur as symmetry-related pairs transporting in opposite 
directions. Moreover, if phase space cannot be decom- 
posed into invariant subsets, e.g., for an ergodic system, 
there cannot be transport even with all symmetries bro- 
ken. 

Up to now we have only considered transport of invari- 
ant sets of the unit cell. For an arbitrary initial distri- 
bution transport is determined by projection onto these 
invariant sets (JO]. Therefore, the location of an initial 
distribution within an invariant set is irrelevant. This 
applies also to the location within the temporal unit cell, 
i.e., to the question of phase dependence discussed in |l2| . 
In particular, in case that the plane p = 0, < x,t < 1 
is completely within the chaotic invariant set, any initial 
condition restricted to this plane will result in the same 
average transport. We now understand how a Hamilto- 
nian ratchet makes particles initially at rest (p = 0) move 
with a predetermined mean velocity as, e.g., in Fig. la. 

We have checked Eq. (|^) numerically for a continu- 
ously driven system || . We determined the areas A{ and 
winding numbers Wi for the regular islands shown in the 
Poincare section of Fig. |l|b as well as (T a ) and (T&) for 
the limiting KAM-tori, yielding v ch = 0.092 ±0.011. The 
error estimate includes the uncertainty in the location of 
the bounding KAM-tori and the contribution from ne- 
glected small islands. The result is in agreement with 
the value i> c h = 0.082±0.002 determined with much more 
computational effort from the spatial distribution of 10 4 
trajectories, started with p = (Fig. |l]a). 

As a minimal model for directed chaotic transport in 
Hamiltonian ratchets, we propose a kicked Hamiltonian 
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H{x,p,t)=T(p) + V(x)^26(t-n). (4) 

n 

It reduces the dynamics to a map for position and mo- 
mentum X n+ i = x n + T'{p n ), p n+1 = p„ - V'(x n+ i), 
just after the kick. As an example we take a symmetric 
potential V(x) = (imod 1 — l/2) 2 /2 and an asymmetric 
kinetic energy T(p) = \p\ + 3 sin(27rp)/(47r 2 ). We con- 
sider the dynamics on a cylinder with transport along 
the x-axis and p G [—1/2, +1/2) being a periodic vari- 
able. Fig. |]d shows the Poincare section for one unit 
cell. There are only two major invariant sets — a chaotic 
sea and a regular island centered around a periodic orbit 
with winding number w rcg = — 1. According to Eq. (|l]), 
transport of the full phase space vanishes identically be- 
cause of the periodic momentum variable. Applying the 
sum rule (|^) the contributions to transport from the two 
invariant sets cancel exactly, 

A ch Wch + Acg W rcg = . (5) 

We find the transport velocity of the chaotic component 
as v ch = / reg /(l - /reg), where / rog = A Ieg /(A Ieg + A ch ) 
denotes the relative area of the regular island. From 
Fig. @d, / reg = 0.117 ± 0.001, thus v ch = 0.133 ± 0.001 
in agreement with u cri = 0.1344 ± 0.0003 from the spatial 
distribution of Fig. |l|c. 



300 





class. , - - 




^-^lCl/256 




- — Ti=1/64 




h^1/8 



t 2000 



FIG. 2. Mean position vs time for a wavepacket in system 
(^) initialized as the momentum eigenstate with p = for 
various values of h (full lines). For decreasing h the classical 
prediction v c h t (dashed line) is approached. 

In order to extend our concept of directed trans- 
port in Hamiltonian systems to quantum ratchets, we 
first demonstrate by a numerical example that quantum 
Hamiltonian ratchets can work. Fig. || shows that the av- 
erage velocity of a wavepacket initialized in the chaotic 
sea varies between for large and the classical value u crl 
for small values of h~. We explain this behavior in the 
following. 

In analogy with our approach to classical transport, 
we consider the invariants of the quantum dynamics, the 
stationary states of the time-evolution operator over one 
period, i.e., U = e - lV / n e - lT / h f or the kicked Hamiltonian 
Eq. (§. They satisfy U\(/) a .k) = exp(-27rie Q (fc)) \<p a , k ), 
with the quasienergy e a (k) £ [0,1) Similarly, spatial 
periodicity implies (j) a ,k{x + l,t) — exp(27rifc) 4> a ,k{x, t) 
with quasimomentum k 6 [0, 1) where h is chosen ra- 
tional for systems periodic in p. Quantum transport is 



related to the expectation values in the stationary states 
v a ,k = {{<f>a,k\v\4>a,k)) of the velocity operator v = T'(p), 
where ((•)) = J Q dx f Q dt(-). Using a generalization of 
the Hellman-Feynman theorem to time-periodic systems 
Jl3| , we express velocities by band slopes as 

v a ,k = de a (k)/dk . (6) 

This allows to discuss quantum transport in terms of 
spectral properties. Examples for quasienergy band spec- 
tra are shown in Fig. ^| together with the correspond- 
ing velocity distributions. The semiclassical regime is 
characterized by the existence of two different types of 
bands and corresponding eigenstates Jl4] , p^[ : Bands per- 
taining to regular states appear as straight lines in the 
spectrum, while the chaotic bands show oscillations and 
wide avoided crossings among themselves. Associating 
the terms chaotic and regular with the bands is supported 
by the Husimi representations of the corresponding eigen- 
functions (insets in Fig. ||a) . The new aspect introduced 
into this picture by directed chaotic transport is the over- 
all slope of the chaotic bands. 
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FIG. 3. Distribution of quantum velocities v obtained ac- 
cording to Eq. (|) for (a) h = 1/128 and (b) h = 1/4. The 
right insets show band spectra for h = 1/32 and h — 1/4, 
respectively. The smaller insets are Husimi representations 
of characteristic wave functions together with the border of 
the classical regular island. In (a) the regular and the chaotic 
wave functions can be associated with the two peaks of P(v) 
centered around u rcg = — 1 and v c h. 

Only on a coarse quasienergy scale, the two sets of 
bands appear to cross. On a sufficiently fine scale, all 
crossings are avoided. Consequently the actual bands 
change their character between regular and chaotic at 
each of the narrow crossings and have no overall slope. 
Switching from the latter ("adiabatic") to the former 
( "diabatic" ) viewpoint is a well-controlled procedure [^4| . 
Formally, this behavior of the bands can be described in 
terms of their winding number (average slope) with re- 
spect to the periodic (e, /c)-space: In the adiabatic as 
well as in the diabatic case, all quasienergy bands must 
close after an integer number of periods in the e and k 
directions, so that their winding number w must be ratio- 
nal. Clearly, in the adiabatic case, all winding numbers 
are zero. Going from the adiabatic to the diabatic case 
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amounts to a mere reconnection of bands at the cross- 
ings, preserving the sum of winding numbers. Thus it 
must be zero also in the diabatic representation, 



E-i res) 



0. 



(7) 



This is the quantum-mechanical analogue of the classical 
sum rule M), Because of the localization of the regu- 
lar states on tori inside the regular island, the winding 
number of the regular bands in (e, /c)-space is in the semi- 
classical limit identical to the winding number in (x, t)- 
space of the central periodic orbit, i.e., wi r ° E ' = w lcg . 
Moreover, in this limit, the fractions of regular and 
chaotic bands correspond to the relative phase-space vol- 
umes / reg and 1 — / reg , respectively. We therefore ob- 
tain from Eq. (]?]) the mean slope of the chaotic bands, 
yj(ch) _ j rcg /(i _ y rcg ) = Vchj as the classical drift veloc- 
ity. This is confirmed in Fig. ||a. 

The asymptotic quantum transport velocity for a given 
initial wavepacket \tp) is an average of band slopes 
weighted with the overlaps \ (ip\(f> a ,k)\ 2 ■ We can now ex- 
plain our observations in Fig. |^: For Ti <C ^4 r egj an initial 
wavepacket prepared in the chaotic region of a single unit 
cell of the extended system is a superposition of chaotic 
eigenfunctions from the entire band spectrum. Conse- 
quently, its drift velocity is given by the mean slope of 
the chaotic bands and thus by the classical value v c h- In 
contrast, for h ^> A rcg , there are no states restricted to 
the regular or the chaotic set and hence quantum trans- 
port does not correspond to classical transport in this 
regime. 
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FIG. 4. Form factor K{x,t) for h = 1/128 (thick line) 
at (a) t = 30 and (b) t = 111, the Heisenberg time of the 
chaotic component. The two distinct peaks centered around 
Dreg and v c h are the fingerprints of directed transport. The 
time-dependent width of the chaotic peak follows the classical 
distribution P' ch '(x,t) of chaotic trajectories (thin line). 

Our analysis based on the winding numbers can be ap- 
plied to predict the mean quantum transport in the semi- 
classical regime from the classical value. The band spec- 
tra, however, encode more detailed information about 
quantum transport. It can be extracted by a double 
Fourier transform e — > t, k — > x (discrete position in 
units of the spatial period) and subsequent squaring of 
the spectral density, translating two-point correlations 
in the bands into the entire time evolution of the spa- 
tial distribution on the scale of the spatial period. A 



formal definition of the resulting generalized form factor 
K(x, t) and further details are found in JlGj . A semiclas- 
sical theory for the form factor in Hamiltonian ratchets, 
which will be published elsewhere, requires to account 
for the simultaneous presence of regular and chaotic re- 
gions in phase space. It relates the form factor K(x, t) to 
the respective contributions of these invariant sets to the 
classical spatio-temporal distribution P(x,t) (Fig. 4). 
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